Finite temperature crossovers and the quantum Widom line near the Mott transition 
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We explore the recently established quantum critical region above the critical end-point in the 
half-filled Hubbard model. The resistivity displays characteristic quantum critical scaling near the 
instability line which we interpret as the quantum Widom line, in a close analogy with the Widom 
line in classical phase transitions. This line is determined by the minimum curvature of the free 
energy functional and represents a continuation of the line of the first order phase transition. The 
quantum Widom line, however, does not require the existence of the first order phase transition. 
It is rather a consequence of the zero temperature quantum phase transition, which in our case is 
masked by the phase coexistence region. Along the quantum Widom line the density of states is 
reduced, but the insulating gap is not yet fully open. We contrast this line with the Fermi liquid 
coherence line which is the boundary of the existence of long-lived quasiparticles. 



I. INTRODUCTION 

Strongly correlated materials exhibit a variety of 
phases whose properties often lack a complete micro- 
scopic understanding.^ While transport and thermody- 
namic properties near the quantum critical point depend 
on details of the system under consideration, at higher 
temperatures they share many common features. A num- 
ber of scaling hypotheses, based on the experimental data 
or theoretical models, are proposed to describe the "high- 
temperature" region of the phase diagram dominated by 
the incoherent excitations .^^^ 

Recently we have established a remarkable scaling 
of the resistivity curves above the critical end-point in 
the frustrated Hubbard model around the instability 
line^ which is determined by the minimum curvature of 
the free energy functional.^ The resistivity around this 
line exhibits a characteristic "fan-shaped" form, surpris- 
ingly similar to some experimental findings,^ reflecting a 
gradual crossover from metallic to insulating transport. 
The scaling behavior in this high-temperature crossover 
regime is interpreted as an evidence of a hidden Mott 
quantum criticality, and it is argued that it incorporates 
universal features of transport in various materials. 

In this paper we explore different regimes in the dy- 
namical mean field theory (DMFT) phase diagram of 
the single band Hubbard model in the paramagnetic 
phase and carefully characterize relevant crossover lines. 
The instability lin^* separating the metallic and insu- 
lating transport is interpreted as a quantum Widom 
line (QWL). It is analogous to the Widom line at the 
crossover between liquid-like and gas-like behavior in su- 
percritical fluids.^ The quantum Widom line represents 
a continuation of the line of the first order phase transi- 
tion and roughly follows the insulating spinodal on the 
metal-insulator transition (MIT) phase diagram. Along 
the quantum Widom line the density of states is reduced, 
but the insulating gap is not yet fully open. We contrast 



this line with the Fermi liquid coherence line which is 
the boundary of the existence of long-lived quasiparticles 
and follows a direction of the metallic spinodal. 

The remaining part of the paper is organized as fol- 
lows. Section II contains the phase diagram featuring 
various transport regimes and crossover lines. In Sec- 
tion III an explicit procedure for a determination of the 
quantum Widom line is given and its physical meaning 
is further illustrated. Section IV shows contour plots of 
the resistivity and /3-function in the whole phase diagram 
and Section V contains the conclusions. 



II. PHASE DIAGRAM AND CROSSOVER 
LINES 

We consider a single band Hubbard model at half- 
filling 

H = - Uj (cj^Cj^ + c.c.) + "^Uni^riii, (1) 

<i,j>a i 

where cj^ and Qct are the electron creation and annihila- 
tion operators, n^^r = cla^ia^ ^ij is the hopping amplitude, 
and U is the repulsion between two electrons on the same 
site. We use a semicircular density of states, and the cor- 
responding half-bandwidth D = 2t\^ set to be our energy 
unit. We focus on the paramagnetic DMFT solution, 
which is formally exact in the limit of large coordination. 
The DMFT is a unique method which studies the Fermi 
liquid low temperature region, the high-temperature in- 
coherent transport and the Mott MIT within a single 
theoretical framework, and it is most reliable precisely 
at h igh temperatures when the correlations are more 
local.^^^^ To solve the DMFT equations we use the it- 
erated perturbation theor}^^^ (IPT) and cross-check our 
results with numericall y exac t continuous time quantum 
Monte Carlo (CTQMC)"^^, as shown in the Appendix. 
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The phase diagram in the U — T plane, Fig. [T] shows 
several regimes and crossover lines. The Fermi liquid 
(FL) region is restricted to small interaction U and low 
temperatures. The boundary of this region is determined 
by a temperature up to which the resistivity follows usual 
p oc behavior (pink dotted line). This happens ap- 
proximately until the resistivity reaches a value O.lp^^^^, 
where p^^^^ is the maximal metallic resistivity in the 
semiclassical Boltzmann theory, which corresponds to the 
scattering length of one lattice spacing. This criterion 
defines the boundary of the FL region and the corre- 
sponding temperature Tfl (see also Section IV). Above 
the FL crossover line the excitations gradually become 
more incoherent until the systems enters deep into the 
bad metal region with p > p^^^^ . It is interesting that the 
Drude-like quasiparticles exist for temperatures substan- 
tially larger than Tp^,, up to the resistivity p ~ pMott^^ 
similar to the doped Hubbard model.^I^ In the half-filled 
model the metallic resistivity can reach values an order 
of magnitude larger than the Pmou-^^^ opening of 

the pseudogap at the Fermi level.^l 

It was recognized early orP^^^ that there are two dis- 
tinct crossover lines in the DMFT phase diagram. Apart 
from the Fermi liquid coherence line which follows the lin- 
ear part the metallic spinodal, there is another crossover 
line - extension of the insulating spinodal, which marks 
a crossover from the bad metal to the activated trans- 
port regime .1^^ By analog}U^^ to the familiar liquid- 
gas transition, we determine the precise location of such 
an "instability trajectory" (green dotted line on Fig. [T]) 
by examining the curvature of the corresponding free 
energy functional.'^ This curvature vanishes at Tc and 
is finite and minimal at T > Tc, along this instability 
line. Consequently, our problem is recast as an eigen- 
value analysis of the corresponding free energy functional 
J^[G(iWn)\ for which the DMFT Green's function solu- 
tion GoMFrii'^n) represents a local extremum j-^^ ' ^^T^ 
The resistivity around such an instability line exhibits a 
remarkably scaling with a power law dependence of the 
scaling parameter,^ as expected in presence of quantum 
criticality. The boundary of the quantum scaling region 
is marked by green dashed lines (see Section III and Ap- 
pendix B). 

In a close analogy with the Widom line which marks 
a crossover between a liquid-like and gas-like behavior in 
supercritical fluids we name the line which determines 
the boundary between a metallic and an insulating be- 
havior a quantum Widom line. There is, however, an 
important distinction: the Widom line is a consequence 
of the first order classical phase transition and it repre- 
sents a continuation of the first order phase transition line 
to higher temperatures. Its quantum analog, the QWL, 
does not require the existence of the first order phase 
transition. It is rather a consequence of the zero temper- 
ature quantum phase transition, which is in the case that 
we study, masked by the phase coexistence region. 

We emphasize that the immediate vicinity T ^ Tc 
of the critical end-point has been carefully studied 
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FIG. 1: Phase diagram of the half- filled Hubbard model. 
Color in the background is the actual result for the resistivity: 
blue color - small resistivity, red color - large resistivity. 

theoreticall}ffS^, and even observed in experiments^, re- 
vealing classical Ising scaling (since one has a finite tem- 
perature critical point) of transport in this regime. In our 
study, we examine the crossover behavior at much higher 
temperatures T ^ displaying very different behavior: 
precisely what is expected in presence of quantum crit- 
icality. The existence of several scaling regimes is well 
known in a study of quantum phase transitions.^^ The 
peculiarity of our scaling analysis stems from a defini- 
tion of the QWL motivated from the physics of classical 
liquid-gas phase transition. 



III. QUANTUM WIDOM LINE 

In this Section we present details of the procedure that 
we use to precisely determine the quantum Widom line. 
We further motivate the proposed term for this crossover 
line between a (bad) metal and an insulator by analyzing 
the resistivity and the density of states along the QWL. 
We explicitly show how several other popular definitions 
for the crossover line essentially reduce to the Fermi liq- 
uid coherence line, in a sharp contrast to the QWL. 



A. Eigenvalue analysis of the free energy curvature 

In order to reveal the proper scaling behavior close 
to the critical end-point, one has to follow a set of tra- 
jectories parallel to "zero field" trajectory Uwl{T)^ and 
dU = U — Uwl{T) is expected to play the role of the 
scaling variable corresponding to a symmetry-breaking 
field favoring one of the t wo co mpeting (metal vs. insu- 
lator) phases. By analog}^^^ to the familiar liquid-gas 
transition, we determine the precise location of such an 
"instability trajectory" by finding the minimum of the 
curvature of the corresponding free energy functionaP^. 
The free energy near such an extremum can be written 
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family of resistivity curves displays characteristic quan- 
tum critical scaling of the form 



p{T,6U) = p^,{T)f{T/Tom). 



(5) 



with To{SU) - |J/7|^^, and p^^(T) - T, with the ex- 
ponent zu ^ 0.56. The scaling parameter To displays 
power law scaling with the same exponents for both scal- 
ing branches and falls sharply as U Uwl^ which is 
consistent with quantum critical scenario and presents 
a stringent test of our scaling hypothesis. In Appendix 
C we show the scaled curves for a broad region on the 
phase diagram around the QWL, and identify the scaling 
region corresponding to the dashed lines in Fig. [l] 



FIG. 2: Eigenvalue analysis: values for A are plotted as a 
function of interaction at three different temperatures. The 
positions of the minima of these curves determine the insta- 
bility line Uwl{T) (quantum Widom line) along which the 
curvature of the free energy is minimal. 



as 



:F[G{iWn)] = ^0 + Tf y SG{iuJm)M, 



m,n 



5G{iuj„) + 

(2) 



where 



M„ 



2Tt2 dG{iojm)dG{iojn) 



(3) 



G^Ge 



and SG{iuJn) = G{iuOn) — GoMFT^'i^n)' The curvature 
of the free energy functional is determined by the lowest 
eigenvalue A of the fluctuation matrix M. As explained 
in Appendix A, A can be obtained from the iterative so- 
lution of DMFT equations. The difference of the Green's 
functions in iterations n and n + 1 of the DMFT self- 
consistency loop is given by 

5G(^+^) {iun) - SG^""^ {iujn) = e-^^(5G(^) {iuj^) (4) 

and therefore A determines the rate of convergence of the 
Green function to its solution. 

An example of our calculations is shown on Fig. |2j 
where the eigenvalues at several temperatures are plot- 
ted as a function of interaction U/Uc- The minima of 
these curves define the locus of the instability trajec- 
tory Uwl{T)^ which terminates at the critical end-point 
(Uc^Tc). The free energy curvature analysis can be used 
to pinpoint the position of the critical end-point from the 
high-temperature data. In this way, one does not need to 
study for this purpose the coexistence region below Tc, 
which in some cases may be accompanied by numerical 
difficulties. 

The resistivity scaling procedure is presented in many 
details in Ref. 8, We first divide each resistivity curve 
by the resistivity along the QWL (the "separatrix" ) 
p^^(T) = p{T,SU = 0), and then rescale the tempera- 
ture for each curve with an appropriately chosen param- 
eter To{SU) to collapse the data onto two branches. The 



B. Density of states and resistivity along the QWL 

The quantum Widom line physically represents a 
boundary between a metallic and insulating transport. 
This can be also seen from the density of states (DOS) 
along the QWL and its vicinity. Middle column on Fig. [3] 
shows the DOS along the QWL for three different tem- 
peratures. While the DOS at the Fermi level is strongly 
suppressed, the gap is not yet fully open. Left column 
shows the density of states in the metallic phase: there 
is a clear Drude peak at low temperatures, which grad- 
ually disappears as the bad metal region is reached by 
increasing the temperature. At larger U the system is in 
the insulating phase with fully open Mott gap, featuring 
activated transport. 

It is also interesting to examine the resistivity along the 
QWL. Its values largely exceeds the Mott limit, yet it dis- 
plays metallic-like but non-Fermi liquid-like temperature 
dependence Pvkl(^) ^ Fig.jlJ Such puzzling behavior, 
while inconsistent with any conventional transport mech- 
anism, has been observed in several strong ly correlated 
materials close to the Mott transitioiP^'^ Our results 
thus suggest that it represents a generic feature of Mott 
quantum criticality. 
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FIG. 3: Density of states around the Widom line. Along the 
Widom line (middle column) the DOS is strongly reduced at 
the Fermi level but the Mott gap is not yet fully open. 



4 




0.04 0.08 0.12 0.16 0.2 



T 

FIG. 4: Resistivity along the quantum Widom line. 



C. Fermi liquid coherence line vs. QWL 

In order to highlight the difference between the quan- 
tum Widom line and se veral other crossover lines of- 
ten used in the literature J^^I^ we plot them together on 
Fig. [5] Light blue line (a) corresponds to p = O.lp^^^^ 
and it is roughly the Fermi liquid coherence boundary 
up to which the resistivity follows p ex T^. At the or- 
ange line (b) the imaginary part of the self-energy at 
two lowest Matsubara frequencies becomes flat. Green 
line (c) is determined as the inflection point of the spec- 
tral density at the Fermi level with respect to the tem- 
perature, d'^A{uj = 0)/(iT^ = 0, and dark blue line 
(d) corresponds to the inflection point of the resistivity, 
d? p{ijj = 0)/dT'^ = 0). All of these four lines essentially 
correspond to a crossover between a Fermi liquid and 
bad metallic transport. In a sharp contrast, the quan- 
tum Widom line, red line (e), determines a boundary 
between the metallic and the insulating transport. 




FIG. 5: Various definitions of the crossover lines. Lines (a)- 
(d) (see the text) essentially correspond to the FL coherence 
line are in sharp contrast to the QWL (e). 



IV. RESISTIVITY AND ^-FUNCTION 

In this Section we plot the resistivity and the /3- 
function in the entire U — T plane. While these plots do 
not directly convey any new physical conclusions, they 
are valuable since they do not include any data analy- 
sis and are, therefore, completely unbiased to any phys- 
ical concept. The DMFT expression for the resistivity, 
p = 1/cf{uo ^0), is obtained from-^-*^ 

/+CO n-\-oo ir 

dev\e)D\e) / {- A\e , uj)) , (6) 
-00 J-00 

where A{e,w) = -^ImG{e,w), v{e) = \J{^t^ - e^)/3. 

^^(^) ~ 2^^^^^ ~ noninteracting DOS, and 

/ is the Fermi function. 

The result is shown on the contour plot on Fig. [g] (resis- 
tivity is again given in units of p^ott)' The white stripes 
follow the lines of equal resistivity and approximately 
separate orders of magnitude in the resistivity. In the 
high-temperature region they emanate from the critical 
end-point (/7c, Tc), which is reminiscent to the quantum 
criticality. 

To specify the scaling behavior, one usually computes 
the corresponding /3-function^ ^[g) = d\ng/d\nT^ with 
g = py^^/p being the inverse resistivity scaling function. 
As we have shown in Ref. |8] it displays a nearly linear 
dependence on In^, and is continuous through SU = 
indicating precisely the same form of the scaling function 
on both sides of the transition - another feature exactly of 
the form expected for genuine quantum criticality. This 
functional form is very natural for the insulating trans- 
port, as it is obtained even for simple activated behavior 
p{T) - e-^^/^. The fact that the same functional form 
persists well into the metallic side is a surprise, especially 
since it covers almost an order of magnitude for the re- 
sistivity ratio. Such a behavior has been interpretecP to 
reflect the strong coupling nature of the critical point, 
which presumably is governed by the same physical pro- 
cesses that dominate the insulator. This points to the 
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FIG. 6: Resistivity calculated in the entire U — T plane. The 
white stripes follow the lines of equal resistivity and approxi- 
mately separate orders of magnitude in the resistivity. 
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FIG. 7: The /3-function calculated in the entire U — T plane 
illustrates the different transport regimes. 



fact that our QC behavior has a strong coupling, 
non-perturbative character. 
The ;5-function, 



I3- 



d\og p 
d log T ' 



I.e. 



(7) 



can be also used to illustrate different transport regimes 
on the phase diagram. The white region on Fig. [7| cor- 
responds to the Fermi liquid regime where p ~ T^. At 
higher temperatures (the green region) resistivity growth 
slows down; precisely in the critical point, {3 diverges 
(yellow) and in its near vicinity, a growth faster than 
quadratic is found (red). The pink line corresponds to 
the resistivity maxima where /3 = 0. At the insulat- 
ing side of the transition, the resistivity is exclusively 
decreasing with increasing temperature and j3 is always 
negative (blue). Deep in the Mott insulator one finds 
the activating regime in which resistivity grows exponen- 
tially with decreasing temperature and at T = 0, /3 goes 
to — oo (purple). 



V. CONCLUSIONS 

In summary, we have explored in detail different 
regimes in the DMFT phase diagram of half-filled frus- 
trated Hubbard model, with a careful analysis of relevant 
crossover lines. The line which separates metallic and in- 
sulating transport is characterized as a quantum Widom 
line in a close analogy with the liquid-gas crossover in 
supercritical fluids, and it is defined as the line of the 
smallest curvature of the free energy functional. This 
line starts from the critical end-point (/7c, T^) on the 
phase diagram, where the curvature of the free energy 
functional goes to zero. The quantum Widom line, how- 
ever, does not require the existence of first order phase 
transition. It rather emanates from the zero tempera- 
ture metal-insulator transition, which is in the case of 



the Hubbard model masked by the phase coexistence re- 
gion. 

We have established a remarkable resistivity scaling 
around the QWL with a power law scaling parameter, as 
expected from the quantum criticality. We believe that 
the proposed concept is relevant for many systems near 
the Mott metal-insulator transition, and it can be tested 
on e.g. various organic Mott systems^ where Tc is suf- 
ficiently bellow room temperature. A stringent test for 
our theory will be on the models where Tc can be signifi- 
cantly reduced^ and where we expect that the quantum 
critical region extends to the lowest temperatures. 
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Appendix A: Eigenvalue analysis 

The minimum curvature of the free energy functional 
for a given temperature is determined as follows. For 
simplicity we concentrate on the Bethe lattice. The 
Ginzburg-Landau free energy functional F{G) in the 
Hilbert space of t he Mat subara Green's functions G = 
G{iuJn) is given by™^ 



F{G) = -Tt^G^^F,mp{G) 



imp 



(G), (Al) 



where the first term is the energy cost of forming the 
Weiss field A = t'^G around a given site, while the second 
term is the free energy of an electron at this site in the 
presence of the Weiss field. 

Close to a local extremum Go, we can Taylor expand 
F{G) in terms of deviation, SG = G — Gq, from this point 

F{G) = F{Go) ^Tt^^6G{icUm)M^^SG{iUn) ^ . . . 



F(Go) +Tt^ SGM6G- 



where 



Mrr 



d^F{G) 



2Tt2 dG{iuJn)dG{iuJm) 



(A2) 



(A3) 



G=Go 
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FIG. 8: Convergence rate in the iterative solution of DMFT 
equations at T/Tc = 2.33 using IPT impurity solver, panel 
(a), and CTQMC impurity solver, panel (b). The dashed 
lines in panel (b) are linear fits to the data. The insets are 
the corresponding eigenvalues determined by the slopes from 
the main panels. 
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FIG. 9: Phase diagram obtained with IPT, panel (a), and 
CTQMC, panel (b). Temperature and interaction are scaled 
by their values at the critical endpoint (T/^^ = 0.046, 
T^^c ^ 03 ^iPT ^ 2.472, U?^^ = 2.3). Red dashed 
line is the instability line U* (T) , full red line is the line of the 
first order MIT, and green and blue dotted lines are left and 
right spinodals. 



We introduce a gradient function 



1 dF{G) 



= MSG 



(A4) 



2Tt^ dG 

and define an iteration-substitution procedure by 

which gives the minimum of the free energy as the itera- 
tion step n ^ 00. In the case of the free energy functional 
given by Eq. (lAll), we have 



^(G) — Gimp{G) — G, 



and in the iterative procedure the Green's function con- 
verges to the minimum of the free energy which is at the 
same time also a self-consistent solution of the DMFT 
equations, given by the relation Gimp{i^n) = G{iuJn)- 

The curvature of the free energy for interaction U and 
temperature T can be obtained as follows. The eigenbasis 
Ga and eigenvalues of matrix M are defined by 



MGa — XqcGcx' 



(A7) 



We can expand 5&^^ as 



(A8) 



where a'a^ are the coefficient s of 8&'^^ in the eigenvalue 
basis. Substituting into Eq. (|A5[), one obtains 



(A6) where 



B„ = -ln(l-A«). 



(A9) 



(AlO) 



For large n the term with lowest = B^q , which cor- 
responds to the lowest eigenvalue Aq,q = A, is dominant 

^qM ^ e-"^-oag^)G«^, n > 1. (All) 

Here ao is the coefficient corresponding to the Green 
function with the lowest eigenvalue A. Now it is obvious 
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that through iterations, the solution G approaches to Go 
exponentiahy along a direction defined by the eigenvec- 
tor of M corresponding to its minimal eigenvalue A. The 
coefficient B^q and the corresponding eigenvalue A are 
then obtained from the slope in the iterative relation 



100 



const — nBao^ (^12) 



which follows from Eq. ( |A9[ ). 

In practice, to obtain A (and thus the curvature of 
free energy), we monitor DMFT loop convergence rate, 
G{iuJo)^'^~^^^ —G{iuJo)''^\ in as many iterations as possible 
and then linearly fit In (G(ia;o)*^"'^"^^ — ^(zcjo)'^"'^) versus 
iteration index n. Here uJq = ttT is the lowest Matsubara 
frequency. For small A, B^q ^ A. We repeat this pro- 
cedure for different values of U at the same temperature 
T to determine Uwl{T) in which X{U)\t is minimal. It 
takes few iterations of the DMF T loop to enter into the 
linear regime given by Eq. ( [aT2| . 

With IPT impurity solver, we can use data from several 
tens of iterations to determine the slope 5^0? Fig. [sja). 
The solution with CTQMC impurity solver has a statisti- 
cal error and the number of iterations is limited before the 
difference \G{iuJo)''^~^^^ — G{iuJo)''^-'\ becomes too small 
and acquires a large relative statistical error. Neverthe- 
less, we were able to determine rather precisely the eigen- 
value A and the interaction U*{T) for which it becomes 
minimal. Fig. [sjb). The "instability line" corresponding 
to the minimum curvature of the free energy is shown in 
Fig.jgfa) (IPT phase diagram), and Fig. [^b) (CTQMC 
phase diagram) . Error bars in Fig. |9|b) are estimates of 
the uncertainty in the position of the instability line. 



Appendix B: Boundaries of the quantum critical 
scaling region 

In this section we explain the procedure we used to 
precisely determine the extent of the quantum critical 
region. In our previous work,^ we defined these bound- 
aries to enclose the region in which the mirror symmetry 
of the two resistivity branches is valid, but the scaling 
itself was found to hold in a broader region. To see at 
what distance from the quatum Widom line the scaling 
formula breaks down and the resisitivity curves do not 
collapse anymore, we calculate the resistivity on a dense 
C/,T grid in a broad region around the QWL and after 
scaling the temperature with 



f = T/To 



2\SU\' 



(Bl) 



we plot the results on Fig. [Toj All the resistivity points 
are found to fall on two very narrow branches when T > 
0.02. Below this value of T the resistivity points start to 
spread and the two branches become increasingly wider. 
We use this estimate to obtain the two quantum critical 
region boundaries (the green dashed lines on Fig.jl]) given 
by 
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FIG. 10: Scaled resistivity in a broad region around the Quan- 
tum Widom line. The scaling formula fails below T/Tq ^ 
0.02. 
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Appendix C: Crossover lines in the CTQMC solution 

All the results presented in this paper obtained with 
the approximate IPT calculations are cross-checked with 
the numerically exact CTQMC method. The agreement 
between the two methods is found to be very good and 
all the features of the QMC results are reproduced within 
the IPT. Here we plot different crossover lines, Fig. [TT] 
and use the notation from Fig. |5] The shape and di- 
rection of all these lines are found to be qualitatively 
the same in both approaches. The line (b) is, however, 
found at relatively higher temperatures in the QMC so- 
lution. Also, the light blue line (a) which corresponds to 
p = O.lp^^^^ and which is roughly the Fermi liquid coher- 
ence boundary up to which the resistivity follows p ex , 
is at lower temperatures than in the IPT. We plot here 
two additional lines: (f) the quasi-particle weight at zero 
temperature defined by 



Z = 



duJn 



and (g) the zero temperature magnetic susceptibility x- 
Both quantities are divided by 10 to fit in the temper- 
ature range of the plot and to be more easily compared 
to the crossover lines. It is evident that the coherence 
temperature is roughly proportional to the quasi-particle 
weight at zero temperature, but with the prefactor 0.1, 
^fl(^) ^ 0-1 ^{U). The quasiparticle weight Z is 
weakly temperature dependent for p < pMott • 
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